Thermal conductance of Andreev interferometers 
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We calculate the thermal conductance G T of diffusive Andreev interferometers, which are hybrid 
loops with one superconducting arm and one normal-metal arm. The presence of the superconductor 
suppresses G T \ however, unlike a conventional superconductor, G T /G% does not vanish as the 
temperature T — * 0, but saturates at a finite value that depends on the resistance of the normal- 
superconducting interfaces, and their distance from the path of the temperature gradient. The 
reduction of G T is determined primarily by the suppression of the density of states in the proximity- 
coupled normal metal along the path of the temperature gradient. G T is also a strongly nonlinear 
function of the thermal current, as found in recent experiments. 

PACS numbers: 74.25.Fy, 74.45. +c, 73.23.-b 
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The thermal transport properties of proximity-coupled 
normal metals have attracted interest recently [l|, 12, L| , 
spurred by recent measurements of the thermopower Sa 
of Andreev interferometers 0, [M @ , which are loops in 
which one arm is fabricated from a superconductor (S), 
and the second arm from a normal metal (N). In these 
experiments, which include some of our own, the ther- 
mopower was found to oscillate periodically as a func- 
tion of applied magnetic flux, with a fundamental period 
corresponding to $o = h/2e. In spite of some recent the- 
oretical attempts |l| , these experimental results have not 
been explained satisfactorily. 

Given this problem, we have more recently turned our 
attention to measurements of the thermal conductance 
G T of proximity-coupled systems |7j , which are more eas- 
ily calculated in the quasiclassical approximation. Bezug- 
lyi and Vinokur |2j recently discussed G T for a NSN sand- 
wich where the size of S is smaller than the coherence 
length. They found a strong suppression of G T with de- 
creasing temperature T and oscillations with the applied 
flux $. However, this behavior was determined primar- 
ily by the modification of the thermal conductance of 
the superconducting component of the sandwich arising 
from proximity to the normal metal, an inverse proxim- 
ity effect. In this Letter, we calculate G T arising from a 
normal-metal proximity effect for two types of Andreev 
interferometers relevant to the experiments, the 'house' 
geometry, and the 'parallelogram' geometry shown in Fig. 
1. Our intent is to determine the effect of sample geom- 
etry and finite NS interface resistance on G T . As be- 
fore, we find that the normalized thermal conductance 
G T /Gjf of these devices decreases as the temperature is 
reduced below the transition temperature of the super- 
conductor Q , but saturates at a finite value that depends 
on the transparency of the NS interface and the distance 
of the NS interfaces from the path of the thermal cur- 
rent. The thermal conductance is highly nonlinear, vary- 
ing with the heat current I T approximately as \[7^ for 
intermediate values of I T . Finally, G T oscillates period- 
ically with an applied magnetic flux with a fundamental 



period of h/2e and an amplitude corresponding to a few 
percent of Gjf, with the amplitude of the oscillations in 
the 'house' geometry being in general larger than those in 
the 'parallelogram' geometry. For both the 'house' and 
the 'parallelogram' geometries, the oscillations in G T are 
symmetric in the applied flux. The suppression in G T is 
directly related to the suppression in the density of states 
in the path of the thermal current. 

Our starting point is the kinetic equations for the non- 
equilibrium distribution functions in the diffusive limit 



equilibrium 



d R [ M oo(d n h L ) + Qh T + M 3Q (d^h T )] = 0, (la) 
9 r [M 33 (d n h T ) + Qh L + M 03 (d R h L )] = 0. (lb) 



Here h l and hr are the components of the full matrix dis- 
tribution function h — h^r + Htt 3 , where r° is the iden- 
tity matrix, and r 3 has the form of the Pauli spin matrix 
o~ z . In the normal-metal and superconducting reservoirs, 
iiL and hx have the equilibrium forms 



h L ,T = -z 



, , E + eV\ , , f E-eV 
tanh — ± tanh 



2k B T 



2k B T 



(2) 



where V is the voltage applied to the reservoir. Mq 
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FIG. 1: Schematic of Andreev interferometers with two dif- 
ferent geometries: (a) 'house' and (b) 'parallelogram'. 
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Ms3, M 03 and M 30 are given by 

m oq = ^ [1 + cosh 6> cosh (9* - sinh sinh 0* cosh(2Q(x))] , 
M 33 = i [1 + cosh 6> cosh 0* + sinh sinh 0* cosh(23(x))] , 
M 03 = i sinh sinh 0* sinh(23(x))), 
M 30 = — sinh 9 sinh 0* sinh(23 i (x)). (3) 

where and x are the solutions of the Usadel equations 
Dd A j a (E,R)=0, (4a) 
Dd\9 - — sinh 29 {d n x? + sinh = 0. (4b) 

In writing these equations, we have assumed that there is 
no superconducting gap in the proximity-coupled normal 
metal, and that the phase x is real in the S reservoirs. D 
is the electronic diffusion coefficient, and j s is the spectral 
supercurrent given by 

j s (E, R) = sinh 2 9(E, R)d n x(E, R), (5) 

and Q in Eqn. (1) is related to j s by Q{E,R) = 
-S(j s (E,R)). 

To close the set of equations, we must specify the 
boundary conditions at the NS interfaces and at the 
nodes where multiple one-dimensional normal-metal seg- 
ments meet. In the limit of low transparent NS inter- 
faces, the boundary conditions can be written in the form 

rsinh6»i(<9£Xi) = sir±h0 2 sin(x2 - Xi), (6a) 
r \dg0\ + i sinh 9\ cosh 0i (e^xi )] = 
cosh 6»i sinh 9 2 e l[x2 ~ Xl) - sinh 9 X cosh 2 . (6b) 

where the subscript 1 (2) refers to the normal metal (su- 
perconductor), r = Rb/RN is a parameter that char- 
acterizes the transparency of the interface, where Rb 
is the resistance of the barrier, and Rn the resistance 
of the normal metal per unit length. In the limit of 
r = 0, these boundary conditions reduce to the conti- 
nuity of the 9 function and the phase x across the in- 
terface. At a superconducting reservoir, 9 is given by 
cosh0 so = E/^/E 2 - |A| 2 , where A is the gap in the 
superconductor. At a normal reservoir, 6Vo = 0. The 
major contribution to the phase x m our simulations is 
from the applied magnetic flux. Since the critical current 
in the superconductor is much larger than that in the 
proximity-coupled normal metal, we assume that all the 
phase change occurs across the normal wires, this change 
being applied symmetrically between the two NS inter- 
faces. For example, for an applied flux of $o/2 through 
the area of the Andreev interferometer loop, we set the 
phases x at the two NS interfaces in the Andreev inter- 
ferometers to be — it/2 and ir/2. For a normal reservoir, 



the absolute value of x is meaningless, but the gradient 
d x x must be 0. 

For a node, assuming one-dimensional wires of equal 
cross-section, the boundary conditions reduce to the con- 
dition that the parameter 9 must be continuous at the 
node, and ^ d x 9 = 0, where the sum is over all wires 
emanating from the node. Similar boundary conditions 
apply to the phase x- 

In terms of these parameters, the thermal current is 
given by 

I T = N D J dE E[M o(d n h L ) + Qh T + M 30 (d M h T )], 

(7) 

where No is the density of states at the Fermi energy. 
It should be noted that the integrand in the equation 
above is simply E times the term in the brackets in Eqn. 
(la), so that Eqn. (la) is simply a statement of the 
conservation of spectral thermal current. 

To determine the thermal current for the Andreev in- 
terferometers, one must solve the coupled systems of dif- 
ferential equations (1) and (4) numerically. In general, 
there is a finite supercurrent flowing between the NS in- 
terfaces in an Andreev interferometer in the presence of 
a magnetic field, so that the terms involving the super- 
current in the equations above cannot be ignored. For 
the geometries shown in Fig. 1, however, some simplifi- 
cations can be made. We shall consider first the 'house' 
geometry (Fig. 1(a)), applying a temperature differential 
AT as shown in the figure, so that one normal reservoir 
is at a temperature T while the other reservoir is at a 
temperature T + AT. Since no supercurrent flows in the 
segments connected to the N reservoirs, the phase x in 
these segments is constant. Furthermore, applying the 
boundary conditions for the phase x as discussed above, 
X = at the central node by symmetry (this is also con- 
firmed by the numerical simulations) , so that we can take 
X = in the segments connected to the normal reservoirs. 
Under these conditions, the last two terms in Eqn. (la) 
vanish. The first integral of the resulting equation is a 
constant, which we denote K\. Integrating a second time 
from the left N reservoir (x = 0) to the right N reservoir 
(x = 2L), we obtain 

h L (x = 2L)-h L (x = Q)=K 1 —dx (8) 

Jo ^00 

The left hand side of the equation above is determined 
by the boundary conditions on the distribution function, 
which are given by Eqn. (2) with V = and the temper- 
atures as defined by Fig. 1. The definite integral on the 
right hand side can be evaluated numerically by solving 
the Usadel equations (4). For a small temperature differ- 
ence AT, in the linear response regime, the left hand side 
can be expanded in a Taylor's series to give the thermal 
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conductance G T = I T /AT 



G = - 



N D 
2k R T 2 
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cosh 2 {E/2k B T) 
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(9) 

Figure 2(a) shows G T of the 'house' interferometer as a 
function of the normalized temperature T/T c , for a sam- 
ple with a perfectly transparent barrier (r = 0). The 
length L of one normal-metal segment attached to a nor- 
mal reservoir defines the correlation energy E c = TiD / L 2 
(see Fig. 1). The transition temperature T c of the super- 
conductor is chosen to be 17.18 E c /kB, corresponding 
closely to the values in the experiments. The different 
curves in Fig. 2(a) correspond to different lengths L' of 
the normal-metal segments connected to the supercon- 
ducting reservoirs. G T differs appreciably from its nor- 
mal state value only at temperatures below E c /kB, 
and saturates as T — > 0. The overall decrease in G T /Gjj 
depends on the length of the segment L', the decrease 
being the greatest for the smallest value of L'. It should 
be noted, however, that even for the smallest value of L', 
G T /Gjv does not vanish at the lowest temperatures, in 
contrast to the case of a pure superconductor ^1 • 

It is well known that the electric conductance G of 
a proximity-coupled normal metal shows the so-called 
reentrance behavior |l3j : as one decreases the temper- 
ature below T c , G first increases and then decreases, ap- 
proaching its normal state value Gm as T — > 0. This 
non-monotonic temperature dependence is a result of the 
competition between two effects. One is the penetration 
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FIG. 2: G T of the 'house' interferometer ((a) and (b)) and the 
'parallelogram' interferometer ((c) and (d)) as a function of 
the normalized temperature T/T c . (a) and (c): for different 
L/L'; (b) and (d): for different r. The insets in (a) and (c) 
show the oscillations of Gt as a function of the phase <f>. 



of superconducting correlations from the S side into the 
N side, which increases the value of G; the other is the 
reduction of the density of states of electrons, which re- 
duces the value of G. The reduction in G T , however, is 
directly associated only with the suppression of the den- 
sity of states N(E, x) in the proximity-coupled normal 
metal 01- This leads to a reduction in the number of 
quasiparticles that can carry a thermal current. Fig. 3(a) 
shows the normalized N(E, x) = cosh(!ft(6>)) cos(3(i9)) for 
the 'house' interferometer along the path of the thermal 
current from x = to x = 2L, for r = and L/L' = 10. 
The suppression is maximum at the node (x = L); how- 
ever, unlike a superconductor, N{E,x) never vanishes 
below the gap. Hence, G T jG T N is suppressed, but re- 
mains finite as T — > 0. In a normal metal, G^ and Gjv 
obey the Wiedemann- Franz law, Gjy oc GatT In 
a proximity-coupled normal-metal, as in a superconduc- 
tor, the Wiedemann-Franz law is violated. Physically, 
this violation results from the enhanced conductance as- 
sociated with the superconducting correlations induced 
in the normal metal, which do not carry any thermal 
current. 

Any parameter that affects N(E,x) will modify G T . 
One parameter is clearly L', the distance between the S 
reservoir and the path of the thermal current, as shown 
in Fig. 2(a). A second parameter is the interface resis- 
tance r. A larger value of r would reduce the proximity 
effect in the normal metal, leading to a smaller change 
in G T . Figure 2(b) shows G T as a function of T/T c for 
the case L/L' — 10, for a number of different values of 
r. As expected, the smaller the value of r, the greater 
the reduction in G T . As both L' and r affect G T , under- 
standing the dependence of G T on L' and r is crucial for 
a quantitative comparison to the experimental results. 

The third parameter that affects N(E, x) is the applied 
magnetic flux. Figure 3(a) shows N(E,x) for zero phase 
difference = between the two S reservoirs; Fig. 3(b) 
shows N(E,x) for <\> — 7T, corresponding to a flux $o/2 
through the interferometer loop. For $o/2, N(E, x) is es- 
sentially flat along the entire normal-metal segment at a 
value corresponding to the normal-metal density of states 
A^q. Consequently, G T regains its normal state value G^ 




FIG. 3: Density of states N(E, x) for the 'house' interfer- 
ometer along the path of the thermal current for r = and 
L/L' = 10. (a) 4> = 0; (b) 4> = jr. 
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FIG. 4: Solid lines represent the strongly nonlinear G T as a 
function of the thermal current I T for the 'house' interferom- 
eter. Dashed lines illustrate the dependence of G T at 
intermediate values of the thermal current. 

at <j> = 7r, so that G T shows full scale oscillations as a 
function of magnetic flux. This can be seen in the in- 
set of Fig. 2(a), which shows G T as a function of <f> for 
r = and L/L' = 5 at a temperature of T/T c = 0.029. 
The observation of full-scale oscillations is directly re- 
lated to the symmetry of the 'house' interferometer. We 
note that the oscillations are symmetric with respect to 
magnetic field. If the central node were not exactly mid- 
way between the two S reservoirs, there would still be a 
suppression of N(E, x) at the node even for </) = %, and 
hence G T would never approach Gjj for any value of (j>. 
Hence the amplitude of the oscillations are reduced if the 
geometry of the 'house' interferometer is not symmetric. 

We now turn to the 'parallelogram' interferometer, 
shown in Fig. 1(b). In this case, in addition to the two 
normal-metal segments of length L connected to the N 
reservoirs, we also have a normal-metal segment of length 
L" in the path of the thermal current. We use a typical 
experimental value of L" /L = 0.66 for all the simula- 
tions in this Letter. Since this segment lies between the 
S reservoirs, a supercurrent may flow in this segment, so 
that the last two terms in the square brackets in Eqn. 
(la) cannot in general be ignored. However, since V in 
the S reservoirs is 0, hx — at the S reservoirs, so that 
the terms involving hx are small and can be ignored. 
This assumption is also supported by detailed numeri- 
cal simulations Q. Under these conditions, the equation 
for G T is again given by an expression similar to Eqn. 
(9), except that the integral of 1/Mqo is over the three 
normal-metal segments along which the thermal current 
flows. 

Figures 2(c) and 2(d) show G T for the 'parallelogram' 
interferometer as a function of L/L' and r respectively. 
The dependence on L/L' and r is qualitatively similar 
to that of the 'house' interferometer. However, unlike 
the case for the 'house' interferometer, the oscillations of 
G T as a function of the phase 4> (shown in the inset of 
Fig. 2(c)) do not reach their full scale amplitude, i.e., G T 
is less than G^ for all values of cj), because there would 



always be a suppression in the density of states around 
the two nodes for this geometry. We also note that the 
oscillations are symmetric with respect to <p (and hence, 
with respect to an applied flux); if the last two terms 
in Eqn. (la) were not negligible, we would also have 
small antisymmetric contributions to G T , since both Q 
and M30 are antisymmetric in the applied flux. 

In experiments on the thermal conductance of Andreev 
interferometers, G T was found to be a strongly nonlinear 
function of the thermal current I T 0, • In order to 
investigate the dependence of G T on I T , we do not take 
the linear response limit of Eqn. (8), but numerically 
calculate K\ for specific values of AT across the inter- 
ferometer. G T is still defined by the ratio I T /AT, but 
it is no longer given by the linear response result Eqn. 
(9). Figure 4 shows the calculated G T as a function of 
I T for the 'house' interferometer (the 'parallelogram' in- 
terferometer shows similar behavior). If G T were linear 
for small I T , we would expect a curve which had zero 
slope at I T — 0, which is clearly not seen in the fig- 
ure. At intermediate values of the thermal current, G T 
varies as VT^, shown by the dashed line in the figure. 
This nonlinear dependence of G T is also observed in the 
experiments and presents a problem in experimentally 
defining the linear-response thermal conductance of the 
sample Q. 
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